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Abstract. This is a theoretical investigation on the formation of the linearly polarized line spectrum of ionized 
cerium in the sun. We calculate the scattering line polarization pattern emergent from a plane-parallel layer 
of Ce II atoms illuminated from below by the photospheric radiation held, taking into account the differential 
pumping induced in the various magnetic sublevels by the anisotropic radiation held. We find that the line 
polarization pattern calculated with this simple model is in good qualitative agreement with reported observations. 
Interestingly, the agreement improves when some amount of atomic level depolarization is considered. We find 
that the best fit to the observations corresponds to the situation where the ground and metastable levels are 
depolarized to about one fifth of the corresponding value obtained in the absence of any depolarizing mechanism. 
One possibility to have this situation is that the depolarizing rate value of elastic collisions is exactly D = 10 6 s _1 , 
which is rather unlikely. Therefore, we interpret that fact as due to the presence of a turbulent magnetic field in 
the limit of saturated Hanle effect for the lower-levels. For this turbulent magnetic field we obtain a lower limit 
of 0.8 gauss and an upper limit of 200-300 gauss. 
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1. Introduction 

When observed towards the limb, the solar spectrum is 
linearly polarized and it shows a surprisingly complex po- 
larization structure, quite different from what might be 
intuitively expected from the standard Fraunhofer spec- 
trum (Stenflo & Keller 1996, 1997). 

One of the striking characteristics of this so-called sec- 
ond solar spectrum is the preeminence of polarization sig- 
nals from several minority species. For example, Ti I and 
Cr I have several multiplets and isolated lines that show 
relatively high fractional polarization values; the same is 
true for several, innocent looking, lines of rare earth ele- 
ments that are hardly noticeable in the Fraunhofer spec- 
trum. This paper deals with the formation of the very 
interesting second solar spectrum of a rare earth: Ce II. 

Rare earths have configurations with more than one 
open shell that yield a rich and complex level structure. 
Ce II is an extreme example; its two lower configurations 
AfhdQs and 4/5<i 2 produce 46 levels within 6000 cm -1 
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from the ground level (Martin, Zalubas, & Hagan 1978), 
and 96 levels within 10000 cm -1 . Such a wealth of low- 
lying levels gives rise to a large number of resonant or 
quasi-resonant transitions, and the population of cerium 
spreads among many metastable levels. This fact, and the 
relatively low abundance of the element (cvce = 1-58 in the 
usual logarithmic scale with an = 12; Asplund, Grevesse 
& Sauval 2005), makes resonant lines of cerium form deep 
in the atmosphere. 

An estimate of the height of formation H at heliocen- 
tric angle 9 of the center of a typical resonance line of Ce 
II (the main ionization stage for solar atmospheric condi- 
tions), is obtained from the expression 
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where we use standard notation and A/"o is the number den- 
sity of hydrogen atoms at the reference height z = 0. For 
ionized cerium, the degeneracy of a typical excited level 
radiatively connected to the ground level is g u — 8 and 
the partition function at the effective temperature of the 
sun is u(T = 5800K) « 230 (van Diest 1980; Irwin 1981). 
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Therefore, a typical visible line (say, at wavelegth A 
5000 A , with Doppler width Av D « 1.65 x 10 9 s" 1 and 
Einstein coefficient for spontaneous decay A u \ ss 10 s s _1 ), 
in an exponentially stratified atmosphere with a scale 
height H w 130 km and 7V = 1.2 x 10 17 cm" 3 (Vernazza, 
Avrett, & Loeser 1976), forms at H w Wlog(7//x) km. 
This is H w 250 km at disc center, and _ff «550 km near 
the limb at = 0.1, i.e., 270 km below and 30 km above 
the photospheric minimum of temperature of the classical 
one-dimensional semi-empirical models, respectively. 

2. Scattering Line Polarization in Ce n 

An elementary analysis of scattering line polarization in 
Ce II may be done along the lines of Manso Sainz & Landi 
Degl'Innocenti (2002, 2003; Papers la, b from now on) for 
Ti i. We consider the excitation state of the atoms in a 
layer located at the top of the solar atmosphere and illu- 
minated from below by the photospheric radiation field. 
Furthermore, we assume that the radiation field is axisym- 
metric with respect to the local vertical and that the polar- 
ization degree is very low (quiet sun case). The excitation 
state of the atoms is then described by the independent 
populations Nm of the (2 J+ 1) magnetic sublevels of each 
atomic level having total angular momentum J, where the 
quantization axis for angular momentum is chosen along 
the symmetry axis 1 . An equivalent and more convenient 
description may be given in terms of the elements (statis- 
tical tensors) 

M=-J V 7 

(2) 

where K — 0, 2, ...,2J — 1 (Ce n levels have half-integer 
total angular momentum). Note that p® is proportional to 
the total population of the level, while p\ is the so-called 
alignment of the level; note also that odd- if elements van- 
ish since no population imbalance between +Mj and —Mj 
sublevels can be created under the hypothesis in this re- 
search. As will become apparent shortly, the emission, ab- 
sorption and dichroism coefficients acquire a very simple 
and transparent form once expressed in terms of the pg 
elements. 

We consider a fairly complete atomic model of Ce n 
which has been compiled from the Database on Rare 
Earths at Mons University (DREAM 2 ; Palmeri et al. 
2000; Zhang et al. 2001). This model has 500 levels and ac- 
counts for 16033 radiative transitions between them. The 
levels have total angular momenta spanning from J = 1/2 
to 15/2, and in all 2038 p$ elements are necessary to de- 
scribe the excitation state of the atom. 

Given an arbitrary level i with total angular momen- 
tum J, the relaxation rate of the statistical tensor p^(i) 

We point out that cerium has no nuclear spin and hence, 
no hyperfine structure. 
2 http : //www. umh. ac .be/- astro/dream. shtml 



due to spontaneous emission from level i towards lower 
levels I is (see Landi Degl'Innocenti & Landolfi 2004) 

-f#{t) y £ l A it , (3) 

i 

while the transfer rate due to spontaneous emission into 
the same level from upper levels u with total angular mo- 
mentum J u is 
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Similarly, the relaxation rate of p$ (i) due to absorptions 
from level i towards upper levels u is 
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and the transfer rate due to absorptions from lower levels 
I to level i, is 

tKt 
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In Eqs. ©-I©, Jo and J5 2 are spherical components of the 
tensor that describes the radiation field illuminating the 
atoms. Neglecting the role of the (typically low) degree 
of polarization of the radiation field on the population of 
atomic sublevels they read: 

4= \ J\'m, (7) 

Jl= -^fd^'W^-lW)], (8) 

where I(p!) is the intensity and p! the cosine of the angle 
that the radiation beam forms with the vertical. 

The balance of Eqs. ©-© establishes the excitation 
state of the atom in statistical equilibrium under a given 
illumination (stimulated emission can be safely neglected 
and has not been considered). The result is an algebraic 
system of equations in the unknowns p$ . This system is 
transformed into a non singular one by substituting one 
of the redundant equations (for example, the equation for 
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Pg of the ground level), with the equation expressing the 
conservation of particles J2 i P(f(i)*/2Ji + 1 = 1. 

The intensity and linear polarization emissivities, 
along a p direction, of a spectral line corresponding to the 
decay from a level u to a level £ may be simply expressed 
as 



line 
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where eo = Qw/4%) A u i^/2J U + 1 N (N being the num- 
ber density of Ce 11 atoms), <f> v is the emission profile 

(2) 

and Wj Je is a numerical factor that depends only on the 
quantum numbers of the levels involved in the transition 
(Table 1 in Landi Degl'Innocenti 1984 gives explicit val- 
ues). Reciprocally, in an optically thick medium absorp- 
tion and dichroism (tjq) coefficients must also be con- 
sidered. They depend on the population and alignment of 
the lower level I 
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where rj = (hvjAm) Bi u ^/2J £ + 1 N. 

As in Papers 1 we use a simplified radiative transfer 
model for the formation of the scattering polarization sig- 
nal. The radiation illuminating the Ce 11 ions is approxi- 
mated by a continuum (spectral details are neglected) , and 
expressions 0-© are evaluated from the observed solar 
center-to-limb variation of the continuum intensity (Cox 
2000). This approximation is entirely justified because of 
the weakness of cerium lines in the intensity solar spec- 
trum. The emergent fractional polarization is then calcu- 
lated in the limit of tangential observation, which gives 



eq/e/ - Vq/vi 

1 - [vq/vi)^q/^Y 



(13) 



where ej and rji, which must be evaluated at the surface of 
the plane-parallel slab, have contributions from both, con- 
tinuum and spectral line. Expression Ijl3|l may be further 
simplified taking into account the low polarization level 
(Vq/Vi < 1- 6 q/ 6 i < !)■ Thus, to first order in eg/ej 
and rjQ/rji, the fractional polarization reduces to (Trujillo 
Bueno 2003) 

9. ~ £9. _ HQ. 

I e/ r\i ' 
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Expression (|14Jl may be further developed taking into 
account that rare earths show very weak lines in the solar 
spectrum. Thus, in the limit r^/rff^ < 1, Eq. 
leads to 
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Fig. 1. Main polarization features in the visible line spec- 
trum of Ce II. Shadowed bars and dashed line represent, 
respectively, the fractional polarization in the spectral 
lines and in the adjacent continuum observed near the 
solar limb (p — 0.1), as reported by Gandorfer (2000, 
2002). The heavy line represents the fractional polariza- 
tion as calculated in this work (p ~ * 0). For representa- 
tion purpouses the calculated fractional polarization has 
been scaled by an arbitrary factor a = 0.08 and added to 
the continuum fractional polarization (see Eq. (|15|) with 
-a). 



^line/^c 



where S lmc and B are the line and continuum source 
functions, respectively. The expression between brackets 
in Eq. (|15JI is easy to evaluate once the p[j and p\ ele- 
ments of the levels involved are known because e^ lnc and 
r]\ lno are given in Eqs. (|§1)- (|12[) above and, assuming that 
the continuum source function is the Planck function, we 
have 
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where we assume a temperature T = 5800 K. The ratio 
^lnuy^cont j n (|15|) depends critically on the details of 
the line formation and cannot be calculated without a true 
radiative transfer treatment. Consistently with our simpli- 
fied transfer model, we assume that this ratio is given by 
some arbitrary value a, identical for all the lines, and we 
take the value that better match the observations. More 
generally, a also accounts for any systematic effect (iden- 
tical for all the lines), not considered in our modeling like, 
for instance, the lower polarization degree for on-disk ob- 
servations with respect to the tangential limit, and the 
broadening of line profiles. 

Figure 1 shows the line core polarization in 26 lines 
of the Ce 11 solar spectrum (shadow bars) as reported by 
Gandorfer (2000, 2002; to be referred to as the Atlas here- 
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after). The lines have been selected as those showing a 
clear polarization peak, either above the continuum or as a 
depolarization feature in the Atlas. We have avoided lines 
with severe blends that do not allow a clear assignment of 
the polarization signal. In particular, we have avoided lines 
having blends with species that are known to have notable 
scattering polarization signals (e.g., Ti I, Cr I, molecules), 
or lines that lie on a highly structured polarization back- 
ground (e.g., on the wings of the H- and K- lines of Ca II). 
We point out, however, that sometimes it is possible to as- 
cribe clearly some Q /I signal to a line despite an intensity 
blend. Those lines have also been included in the study. 

The heavy line in Fig. 1 shows the polarization sig- 
nal in each transition calculated according to the anal- 
ysis above (Eqs. l|S|)- (|15|) '). Our simplified scattering po- 
larization modeling is able to reproduce the sign of most 
of the selected lines and roughly the relative Q/I signals 
among lines. A quantitative evaluation of the agreement 
between the observed polarization pattern and the the- 
oretical estimate may be given through the correlation 
factor r — ^Pili/C^Pi^Qi) 1 ^ 2 > where pi and q t corre- 
spond to the observed and theoretical fractional polariza- 
tion in the i-th. line, referred to the continuum. For the 
case shown in Fig. 1, r = 0.72. 

A better modeling of the observed pattern requires 
to include depolarizing collisions and turbulent magnetic 
fields. A detailed treatment of depolarizing collisions is 
not possible because of the lack of data. Some progress is 
indeed going on on the theoretical side (e.g., Derouich et 
al. 2003, 2004; Kerkeni 2002), but no reliable calculations 
exist (to the authors' knowledge) of depolarizing rates in 
rare earth atoms. Nonetheless, much qualitative under- 
standing on their role may be gained by proceeding as in 
Papers I; we just consider one and the same depolarizing 
collisional rate D for all the statistical tensors with K ^ 
in all the levels and we add a rate of the form —Dpq(i) to 
the statistical equilibrium equation of each p$ (i) element 
with K 0. 

Figure 2 shows the fractional polarization calculated 
assuming a depolarizing rate D = 10 6 s _1 for all the 
atomic levels (heavy line). The agreement with the ob- 
served polarization pattern (shadow bars) is better than 
in the D = case (cf. Fig. 1); quantitatively, this is re- 
flected in a larger correlation factor r — 0.83. However, 
the agreement does not improve further for increasingly 
larger depolarizing rates but, on the contrary, it degrades, 
which puts an upper limit to the depolarizing collisional 
rate at about ~ 10 6 s _1 . 

A depolarizing rate D — 10 6 s~ : leaves the atomic po- 
larization of the lowermost levels (ground and metastable) 
of Ce II at about one fifth of the atomic polarization in 
the absence of depolarizing processes. This is interesting 
because a microturbulent magnetic field may also depo- 
larize an atomic level whenever its corresponding Larmor 
frequency vl is of the order of the relaxation rates of the 
level (i.e., the inverse of its mean life-time r) or larger. But, 
differently from the case of collisions, the depolarization 
saturates when vl ^> 1/r, and the atomic polarization 
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Fig. 2. Same as Fig. 1 with a common depolarizing rate 
D = 10 6 s _1 for all the levels of the Ce II model, and 
a = 0.12. 

of the level remains at a constant value that is one fifth 
that of the non magnetic case (see Appendix). Therefore, 
the results in Fig. 2 are consistent with the presence of a 
turbulent magnetic field in the saturation regime of the 
Hanle effect for the lower levels. The lowermost levels of 
our Ce II model atom have mean life times r «l-8x 10 -6 s; 
the turbulent magnetic field required for reaching the sat- 
uration regime of the lowermost levels of Ce n is hence of 
the order of 0.1-0.8 gauss. 

On the other hand, the magnetic field cannot be so 
strong to saturate also the upper levels of the transitions, 
or the fractional polarization pattern would be exactly as 
in the non magnetic case but reduced to one fifth, and 
therefore, r would degrade again to 0.72. The upper levels 
of the transitions considered in this study have mean life- 
times r w 4-6 xlO -9 s; most of these levels saturate with 
fields B > 200 gauss, all of them when B > 300 gauss. 

3. Conclusions 

Rare earths are atoms with a complex level structure 
whose scattering polarization patterns cannot be under- 
stood using simple two level approximations. Moreover, in 
astrophysical environments (and in the solar atmosphere 
in particular), many transitions are excited simultaneously 
and the population and atomic polarization of a level de- 
pend on the balance between all the absorption and emis- 
sion processes involving the level. 

Here, we have considered the generation and distribu- 
tion of atomic alignment in a multilevel atomic model of 
Ce II when it is illuminated by a solar phostospheric-like 
radiation field. We have seen that the ensuing polariza- 
tion of the scattered light is in good agreement with the 
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observed polarization pattern of Ce II. In particular, our 
model is capable of explaining the strongest polarized sig- 
nals, as well as the sign of the polarization signal observed 
in most of the lines selected in this study. 

With a simple analysis of the role of depolarizing col- 
lisions we have shown that the best agreement with ob- 
servations is achieved at depolarizing rates D = 10 6 s _1 , 
thus establishing an upper limit to these collisions. In this 
regime, the lowermost levels of ionized cerium show an 
atomic polarization that is of the order of one fifth of that 
which would be present in the absence of any depolarizing 
mechanisms. We interpret this fact as an indication of the 
presence of volume-filling microturbulcnt magnetic field of 
strength B > 0.8 gauss (the saturation regime of the Hanle 
effect for the lowermost levels) . Since in the presence of a 
magnetic field strong enough to saturate also the upper 
levels the relative polarization of the lines would match 
observations exactly as in the non magnetic case (which 
is not the best fit), we conclude that the magnetic field 
is not saturating the upper levels and hence the magnetic 
field must be B < 200-300 gauss. 

In summary, we have given a theoretical interpretation 
of the second solar spectrum of Ce II and shown that it is 
a promising diagnostic tool to probe the magnetism and 
thermodynamics of the solar atmosphere. However, we 
cannot proceed much further with the elementary meth- 
ods used in this paper. 

Two are the main limitations of this research that must 
be addressed next. First, the full radiative transfer prob- 
lem must be solved, thus allowing for different heights of 
formation of the lines. Second, the role of the magnetic 
field must be included consistently beyond the analogy 
with depolarizing collisions used here. The results of these 
investigations will be presented in forthcoming papers. 
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Appendix A: Saturation Regime of the Hanle 
Effect from a Microturbulent Field 

When a magnetic field is present in an otherwise axisym- 
metric enviroment, it is not possible, in general, to char- 
acterize the excitation state of the atoms only with the 
populations of the individual magnetic sublevels (or equiv- 
alently, through the statistical tensors p$). Taking the 
quantization axis along the direction of the magnetic field, 
the radiation field is characterized through the Jq com- 
ponents (K = 0,2, Q = -K,...,K), where the Q ^ 
elements are azimuthal averages of the intensity of the ra- 
diation field (e.g., Landi Degl'Innocenti & Landolfi 2004). 
Correspondingly, the atomic excitation is described by the 



larger set of statistical tensors Pq , (K = 0, 2, 4, K max , 
Q = —K, ■■■,K\ K max = 2 J or 2 J — 1 for integer or half- 
integer J, respectively), where the Q ^ elements corre- 
spond to quantum coherences between magnetic sublevels 
(e.g., Fano 1957; Landi Degl'Innocenti & Landolfi 2004). 
If the radiation field satisfies the condition Jq /Jq <C 1 
(for, K , Q ^ 0), which is typically the case in stellar at- 
mospheres, then Pq/po <C 1 (for, K,Q ^ 0), and a per- 
turbative treatment of the statistical equilibrium equa- 
tions can be performed (weak anisotropy approximation; 
see Landi Degl'Innocenti 1984; Landi Degl'Innocenti & 
Landolfi 2004, Sect. 10.13). Then, to first order, the sta- 
tistical equilibrium equations in the magnetic field refer- 
ence system do not couple the pq components with pq, 
components having Q ^ Q' . 

Let V and M. be the reference systems with the quan- 
tization axis directed along the vertical and the magnetic 
field, respectively; let R be the rotation taking V into M, 
and Vqq,(R) (Q,Q' = —K,...,K) the corresponding ro- 
tation matrix (e.g., Brink & Satchler 1968). Let [pq]° be 
the alignment of a level in the reference system V, for 
zero magnetic field. In the reference system M., it turns 
out that [po] m = ^ > oo(^)[/ 9 o]°- This is true even if the mag- 
netic field is not zero, since it does not alter [pq] 111 directly, 
and [po] m is decoupled from other [pq] m components with 
Q 7^ (see above). Moreover, if the Larmor frequency ul 
is much larger than the radiative rates IZi of the level, 
this is the only non zero density matrix element because 
[pq] m — > for components Q ^ 0. Finally, passing from 
M to V, = V^ Q {R)*V^{R)[pl) . Averaging over all 
the directions of the turbulent field, [pq] = <5qo[/°o]°/5- 
Therefore, the alignment of a level in the saturated Hanle 
effect regime (vl TZi), is 1/5 of its 'non magnetic' value. 

Now, suppose there is magnetic field regime such that 
the atomic levels can be classified in two types; type A: 
levels in the saturated regime and type B: levels in the 
regime vl CK,. The alignment of levels in set A is 1/5 of 
the value obtained by solving Eqs. The alignment 

of the levels in set B is obtained solving a restricted system 
of equations under this constraint. Moreover, it can also 
be shown that the alignment of the levels in set B is the 
weighted average (with weights 1/5 and 4/5, respectively), 
of the value obtained by solving Eqs.©-©, and of the 
value obtained by solving the same equations under the 
constraint of complete depolarization of the levels of set 
A. 

The rationale for this problem is the following. Because 
of the large life-time difference between metastable and 
excited levels, the lower- and upper-level Hanle effect 
regimes are usually well separated (ti/t u ~ A u i/ Bi u Jq ~ 
100 for two levels conected by an allowed visible tran- 
sition). Therefore, there is a regime of magnetic fields in 
which metastable levels are completely saturated while ex- 
cited levels remain essentially unaffected by the magnetic 
field. In some astrophysical enviroments the microscopic 
fields created by the motion of nearby ions may be large 
enough to saturate the Hanle effect of metastable levels. 



6 



R. Manso Sainz et al.: Second Solar Spectrum of Ce II 



From a physical point of view, the zero (macroscopic) field 
case then corresponds to the solution of the problem just 
stated. 
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